Search for New Particles Decaying to tt in pp Collisions at y/s = 1.8 TeV 
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We use 106 pb _1 of data collected with the Collider Detector at Fermilab to search for narrow- 
width, vector particles decaying to a top and an anti-top quark. Model independent upper limits on 
the cross section for narrow, vector resonances decaying to ti are presented. At the 95% confidence 
level, we exclude the existence of a leptophobic Z' boson in a model of topcolor-assisted technicolor 
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with mass M z , < 480 GeV/c 2 for natural width Y = 0.012 M z >, and M z , < 780 GeV/c 2 for T = 
0.04 M z , . 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk 
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In this letter, we present a model-independent search 
for narrow, vector resonances decaying to ti. This search 
is sensitive to, for example, a Z' predicted by topcolor- 
assisted technicolor |j],|2j . This model anticipates that the 
explanation of spontaneous electroweak symmetry break- 
ing is related to the observed fermion masses, and that 
the large value of the top quark mass suggests the intro- 
duction of new strong dynamics into the standard model. 
It accounts for the large top quark mass by predicting the 
existence of a residual global symmetry SU(3)xU(l) at 
energies below 1 TeV. The SU(3) results in the genera- 
tion of topgluons which we have searched for previously 
in the bb channel p|- The U(l) gives the Z' we search 
for here. In one model g], the Z' decays exclusively to 
quarks (lcptophobic) resulting in a large cross section for 
ti. 

With the z-axis defined along the proton beam, the 
Collider Detector at Fermilab (CDF) coordinate system 
defines <f> as the azimuthal angle in the transverse plane, 9 
as the polar angle, and pseudorapidity 77 as — In (tan |). 
Tracking chambers, immersed in a 1.4-Tesla solenoidal 
magnetic field, are used for the detection of charged 
particles and the measurement of their momenta. The 
precision track reconstruction of the silicon microstrip 
vertex detector (SVX), located immediately outside the 
beampipe, is used for the detection of displaced sec- 
ondary vertices resulting from &-quark decays. Outside 
the SVX is the vertex time projection chamber (VTX) 
which provides further tracking information for < 
3.25. Both the SVX and VTX are housed within the cen- 
tral tracking chamber (CTC), a wire drift chamber used 
to measure charged particle momenta. Electromagnetic 
and hadronic calorimeters, located beyond the CTC and 
superconducting solenoid, measure energy in segmented 
r)-(j) towers out to \r]\ < 4.2. Drift chambers used for muon 
detection reside outside the calorimetry. A more detailed 
description of the CDF detector can be found elsewhere 

Standard model tt production in pp collisions at a cen- 
ter of mass energy of y/s = 1.8 TeV is dominated by qq 
annihilation, while ~10% is attributable to gluon-gluon 
fusion. Once a pair is produced, each of the top quarks 
is expected to decay almost exclusively to Wb. The 
search presented here focuses on the ti event topology in 
which one W boson decays hadronically while the other 
decays to an electron or muon and its corresponding neu- 
trino. The fragmentation of the 6-quarks, as well as the 
hadronic daughters of the W boson, form jets. Accord- 
ingly, ti candidates in this "lepton + jets" channel are 
characterized by a single lepton, missing transverse en- 
ergy, Fj t @, due to the undetected neutrino, and at least 
four jets. Furthermore, a jet resulting from a 6-quark can 
be identified (or "tagged" ) as such by the reconstruction 
of a secondary vertex from the b hadron decay using the 
SVX, or by using the soft lepton tagging (SLT) algorithm 
to find an additional lepton from a semileptonic b decay 

w 

Like the top quark mass measurement [pi, events in- 



cluded in our measurement of the ti invariant mass spec- 
trum must first contain a lepton candidate in the central 
detector region (|ry| < 1.0). This lepton is required to be 
either an isolated electron with transverse energy (E T ) 
in excess of 20 GeV or an isolated muon with transverse 
momentum (P T ) in excess of 20 GeV/c. Events must also 
include at least 20 GeV of Hl T , attributable to the pres- 
ence of a neutrino, as well as at least four jets with \rj\ < 
2.0 and raw E T > 15 GeV. Raw jet energies are the val- 
ues which result from clustering signals in the calorimeter 
towers before any offline jet corrections are applied. To 
increase the acceptance rate for ti events, the require- 
ments for the fourth jet are relaxed such that the raw E T 
must only be greater than 8 GeV with \rj\ < 2.4 in events 
where at least one of the leading three jets is tagged by 
the SVX or SLT algorithms. All jets in this analysis are 
formed as clusters of calorimeter towers within cones of 
fixed radius AR = \J{A-q) 2 + (A0) 2 = 0.4. In 106 pb" 1 
of data, we observe 83 events which satisfy these require- 
ments. 

This method builds upon the techniques developed for 
the top quark mass measurement [^j by fitting each event 
to the hypothesis of ti production followed by decay in 
the lepton+jets channel: 

pp — ► t i £ 

I -W- b 

I q q' (or l~ v x ) 

"W + b 

►■ l + vi (or q q') 

The four-momenta of these 13 objects fully describe a 
ti event. The three-momenta of the charged lepton and 
four jets are measured directly. To compute the energies 
of these objects, the b and b quark masses are taken to 
be 5 GeV/c 2 , the q and q' masses are taken to be 0.5 
GeV/c 2 , and the charged lepton mass is assigned accord- 
ing to its identification as either an electron or a muon. 
The components of transverse momentum for the recoil- 
ing system, £, are measured directly from extra jets in 
the event and unclustered energy deposits that are not 
included in lepton or jet energies. The transverse mo- 
mentum components of the neutrino are computed by 
requiring that the total E T in the event sums to zero. 
While the neutrino is assumed to be massless, its longitu- 
dinal momentum is a free parameter in the kinematic fit 
in which the qq' and Iv invariant masses are constrained 
to equal the W boson mass. We perform a kinematic fit 
to the production and decay of the ti pair as described 
by the decay chain shown above. This fitting procedure, 
which depends on the minimization of a \ 2 expression f| , 
allows the lepton energy, the jet energies and the unclus- 
tered energy to vary within their respective uncertainties. 
The fitted results for these values determine the t and i 
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four-momentum, from which the ti invariant mass (M t j) 
can be computed. To improve the M t j resolution, we also 
constrain the two Wb invariant masses to 175 GeV/c 2 , in 
agreement with the most recent measurement of the top 
quark mass |u|. We use only the four highest E T jets, 
leading to 12 combinations for assigning jets to the b, b, 
and hadronic W daughters. However, because we mea- 
sure only the transverse component of the total energy, 
thereby determining j! T , a two-fold ambiguity in the lon- 
gitudinal component of the neutrino momentum results 
in 24 combinations. We further require that jets that 
are SVX or SLT-tagged be assigned to 6-quarks, thereby 
reducing the number of combinations. 

Electron energies and muon momenta are measured 
with the calorimeter and tracking chambers, respectively 
jyj. A set of generic jet corrections is applied to the 
energies of all the jets in an event to account for ab- 
solute energy scale calibration, contributions from the 
underlying event and multiple interactions, as well as en- 
ergy losses in cracks between detector components and 
outside the clustering cone. These corrections are deter- 
mined from a combination of Monte Carlo simulations 
and data jl2| . The four leading jets in a ti event undergo 
an additional energy correction that depends on the type 
of parton that they are assumed to be in the fit: a light 
quark, a hadronically decaying b quark, or a b quark that 
decayed semileptonically. These parton-specific correc- 
tions account for (a) the differences in the expected P T 
distributions of jets from ti and the shape which was 
assumed to derive the generic jet corrections mentioned 
above, and (b) the energy losses from semileptonic b and 
c-hadron decays. These corrections were derived from 
a study of ti events generated with the herwig Monte 
Carlo program p3| . 

Using Monte Carlo simulations of signal and back- 
ground events, we explored several event selection cri- 
teria in an attempt to optimize our discovery potential 
. Of the 24 possibilities for each event, we select the 
M t i value which corresponds to the configuration with 
the lowest x 2 ■ To reduce the probability of selecting 
configurations with incorrect parton assignments which 
tend to yield artificially low values of M t t, we refit each 
event after releasing the constraint that the Wb invari- 
ant mass be equal to 175 GeV/c 2 and demand that the 
fit for this particular configuration return a value for the 
top quark mass between 150 GeV/c 2 and 200 GeV/c 2 . 
To further reduce incorrect combinations and to increase 
discovery potential for a new particle decaying to ti, we 
apply a y 2 cut. For narrow width ti resonances, simula- 
tion predicts that the width of the M t i spectrum is ~ 6% 
of the resonance mass for cases in which the correct jet 
configuration is selected. For resonances with a natural 
width r that is significantly less than 6% of the nominal 
mass, the CDF detector resolution will dominate and the 
resonances will all have approximately the same shape 
(shown in the inset of Fig. 1 for a mass of 500 GeV/c 2 ). 
At low M t j, the presence of residual events with incorrect 
parton assignments is evident in this figure. 



The selection criteria described above eliminate an ad- 
ditional 20 events from our data sample and the resulting 
M t i spectrum is shown in Fig. [l], along with the expected 
standard model ti and QCD M^+jets background shapes 
normalized to the data. While the non-tt background is 
dominated by VF+jets events, it also includes contribu- 
tions from multijet bb events with one jet misidentified 
as a lepton, Z+jets events, events with a boson pair, 
and single-top production. However, it has been shown 
that the VECBOS |L5| VF+jets shape alone is sufficient in 
modeling the entire non-tt background spectrum || . For 
this analysis, the expected mm-ti background prediction 
of 31.1 ± 8.5 events is calculated as in Ref. ||, but ac- 
counts for differences in selection criteria. We find that 
the M t i distribution of 63 data events is consistent with 
the hypothesis that the spectrum is comprised of stan- 
dard model ti production and the predicted rate of non-ti 
background events, as shown in Fig. |l|. 

Because we cannot present evidence for a narrow ti 
resonance, we establish upper limits on the production 
cross-section for a new vector particle, X, of mass M x 
decaying and to ti. For natural widths T — 0.012M X 
and T = 0.04A/ X , and for each M x between 400 GeV/c 2 
and 1 TeV/c 2 in increments of 50 GeV/c 2 , we perform a 
binned- likelihood fit of the data. To determine the likeli- 
hood function for a given M x and T, we fit the M t j spec- 
trum from the data to the expected Monte Carlo shapes 
for both the ti and QCD H^+jets background sources as 
well as a resonance signal X — > ti which we model using 
Z' — > ti in pythia Q . 

Our analysis is subject to several sources of systematic 
uncertainty on the expected shape of background and sig- 
nal M t i spectra and/or the signal acceptance rate. Treat- 
ing these two types of systematic effect separately, we 
establish the magnitude of each source through a Monte 
Carlo procedure which quantifies the effect of varying 
the source of uncertainty by one standard deviation. We 
determine the uncertainty contributions due to the jet 
E T scale, initial and final state gluon radiation, and the 
non-tt background spectrum using methods described in 
Ref. || . The uncertainty in the measurements of the top 
quark mass |l(| and total integrated luminosity JTt) are 
included in our study of systematic effects, as well as the 
uncertainty due to the choice of parton distribution func- 
tions (PDF) . The remaining sources of systematic uncer- 
tainty considered are all small and include trigger effi- 
ciency, lepton identification, tracking efficiency, z-vertex 
efficiency, and Monte Carlo statistics. The uncertainties 
resulting from jet E T scale and top quark mass are cor- 
related and we conservatively take this correlation to be 
100%. 

The percent uncertainty in ct x • BR{X — > ti} is listed 
in Table | for each of the systematic sources at several 
different resonance masses. The systematic effect due to 
uncertainty in the top quark mass (M top ) is dominant 
at low M x , whereas the effect due to the uncertainty in 
modeling final state radiation dominates at large M x . 
To ensure that our estimates are conservative, the sys- 
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tematic uncertainty is taken to be a constant number of 
pb below the value of a x ■ BR{X — > tt} corresponding to 
the 95% C.L. limit obtained with statistical uncertainties 
only Jl4| . That constant is the estimate of the systematic 
uncertainty at the 95% C.L. limit. Above the same value 
of er x • BR{X — > tt}, we use a systematic uncertainty that 
rises with cr x ■ BR{X — * tt} at the fixed percent rate listed 
in Table § 

For each resonance mass and width, we convolute the 
statistical likelihood shape with the Gaussian total sys- 
tematic uncertainty and extract the 95% C.L. upper limit 
on cr x • BR{X — > tt} which is listed in Table || and shown 
in Fig. ||. The systematic uncertainties increase the 95% 
C.L. upper limit by 27% for M x = 400 GeV/c 2 , but only 
7% (6%) for Af x = 600 (800) GeV/c 2 because statisti- 
cal uncertainties dominate the likelihood. Also shown in 
Fig. H are the theoretical predictions for cross-section 
times branching ratio for a leptophobic Z' with natu- 
ral width L = QM2M Z , and T = 0MM z , ||. At 95% 
confidence, we exclude the existence of a leptophobic top- 
color Z' with mass M z > < 480 GeV/c 2 for natural width 
L = 0.012 M z >, and mass M z > < 780 GeV/c 2 for L = 
0.04 Mz 1 - For larger widths, detector resolution will no 
longer be the dominant factor in determining the Z 1 sig- 
nal shape, so our limits are no longer applicable. 

In conclusion, after investigating 106 pb -1 of data col- 
lected at CDF, we find no evidence for a tt resonance 
and establish upper limits on cross-section times branch- 
ing ratio for narrow resonances. We have used these lim- 
its to constrain a model of topcolor assisted technicolor. 
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QCD background shape 
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Additional acceptance effects 
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TABLE I. The percent systematic uncertainty in 
<7 X ■ BR{X — * tt} from various sources. 
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TABLE II. The 95% C.L. upper limit on the cross section 
times branching ratio for vector particles decaying to ti, as a 
function of mass, for two natural widths. 
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FIG. 1. The observed M t j spectrum (points) compared to 
the QCD iy+jets background (fine dashes) and the total stan- 
dard model prediction including both QCD W+jets and ti 
production (thick dashes). The ti prediction has been nor- 
malized such that the number of events in the total stan- 
dard model prediction is equal to the number of events in 
the data. The inset shows the expected M t i shape resulting 
from the simulation of a narrow, vector X — ► ti resonance 
(M x = 500 GeV/c 2 , T = 0.012M X ) in the CDF detector. 



CDF 95% C.L. Upper Limits for r = 0.012M X 
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FIG. 2. The 95% C.L. upper limits on ct x • BR{X -> ti} as 
a function of mass (solid and open points) compared to the 
cross section for a leptophobic topcolor Z' (thick solid and 
dashed curves) for two resonance widths (r = 0.012 M z > and 
r = 0.04 M z ,). 
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